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Abstract 

We investigate the low temperature thermodynamics and correlation functions of one- 
dimensional spin-1/2 fermions with strong repulsion in an external magnetic field via the thermo- 
dynamic Bethe ansatz method. The exact thermodynamics of the model in a weak magnetic field 
is derived with the help of Wiener-Hopf techniques. It turns out that the low energy physics can 
be described by spin-charge separated conformal field theories of an effective Tomonaga-Luttinger 
liquid and an antiferromagnetic SU{2) Heisenberg spin chain. However, these two types of confor- 
mally invariant low-lying excitations may break down as excitations take place far away from the 
Fermi points. The long distance asymptotics of the correlation functions and the critical exponents 
for the model in the presence of a magnetic field at zero temperature are derived in detail by solv- 
ing dressed charge equations and by conformal mapping. Furthermore, we calculate the conformal 
dimensions for particular cases of correlation functions. The leading terms of these correlation 
functions are given explicitly for a weak magnetic field H <^ 1 and for a magnetic field close to the 
critical field H — J- He- Our analytical results provide insights into universal thermodynamics and 
criticality in one-dimensional many-body physics. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.IK, 05.30.Fk 
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I. INTRODUCTION 



Since the pioneering work in the 60's, 70's and 80's by McGuire, Yang, Lieb, Sutherland, 
Baxter et ai, and of the St Petersburg and Kyoto schools, the study of integrable models 
has flourished into a major activity. Almost without exception, the energy levels are given 
exactly in terms of the Bethe ansatz (BA) equations, from which physical properties can 
)e calculated. This is a hallmark of integrable models that exhibit Yang-Baxter symmetry 
l| . The knowledge and understanding gained from integrable models have greatly enhanced 
progress in the theory of phase transitions and critical phenomena. The most significant 
results achieved to date have been for two-dimensional lattice models and their related 
one- dimensional (ID) quantum spin chains as well as strongly correlated electronic systems 
{2-6]. Integrable models are also known for systems such as Bose- Einstein condensates [?, 8|, 
metallic nanograins [9] and impurity models lOl. 

In general, the BA solution for ID integrable systems is a set of coupled algebraic equa- 
tions. Finding a set of solutions of quasimomenta and spin rapidities {kj,Aj} for the BA 
equations gives the energy V . /cf and the momentum . fc,- of the system. However, the BA 
equations by themselves do not explicitly exhibit any temperature dependence. At zero tem- 
perature, the BA equations in principle give the complete eigenstates of the model. However, 
at finite temperatures, the equilibrium states become degenerate. Thus the thermodynamics 
of these BA solvable models are instead determined by a set of coupled nonlinear integral 
equations called the thermodynamic Bethe ansatz (TBA) equations 12|]. The TBA equations 
are expressed in terms of the dressed energies of different "Fermi seas" that are functions 
of temperature, chemical potential and external magnetic fields. The TBA equations in the 
zero temperature limit, i.e., T — )■ 0, give rise to the so called dressed energy equations which 
describe the band fillings with respect to Zeeman fields and chemical potentials. 

The TBA equations are very difficult to solve in general. They involve an infinite number 
of coupled nonlinear integral equations for spin strings which are quite cumbersome to solve 

rifi nn 

using either analytical or numerical methods |^|6[. Recently, Caux et al [13|, developed 
numerical schemes to solve the TBA equations of the ID two-component spinor Bose gas 
with delta-function interaction. The results obtained by these numerical schemes show an 
insightful interplay between quantum statistics, interactions and temperature in ID inter- 
acting many-body systems. In the context of the quantum transfer matrix method 15| . 
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Kliimper and Patu [16j derived the nonlinear integral equations for the ID Bose and Fermi 
gases with repulsive delta-function interaction. This approach opens up the possibility of 
obtaining the thermodynamics of the continuum models of interacting fermions and bosons 
by taking an appropriate limit for the integrable lattice models. The advance of such ap- 
proaches is the reduction of the infinite number of TBA equations to a finite number of the 
nonlinear integral equations. Where the finite number of the nonlinear integral equations 
for the lattice models can be solved numerically and analytically in certain temperature 
regimes 
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18| . Despite giving high precision numerical thermodynamics, finding the uni- 
versal nature of interacting particles requires further analytical input. Significant universal 
features of ID many-body systems are Tomonaga-Luttinger liquid physics and quantum 
critical phenomena at low temperatures which involve finding essentially universal param- 
eters such as central charges, Luttinger parameters, correlation exponents and dynamical 
critical exponents. All studies of these universal parameters call for mathematical analysis 
and analytical derivation. 

Some progress has been made to derive low temperature analytic results for BA solvable 



models. For example, Mezincescu et al. 19|, |20| obtained the free energy of spin chains at 



low temperatures under a small magnetic field by using the Wiener-Hopf technique. Johnson 
and McCoy [21] obtained the leading temperature dependent terms in a low temperature 
expansion of the free energy for the massive regime of the Heisenberg model. Filyov et al. 



22| gave an exact solution to the s-d exchange model expressed as a series in terms of the 



temperature. Some analytical results for the TBA equations of ID many-body systems are 



restricted to the ground state (T = 0) in the strong coupling limit (c ^ 1) |23l-l25|. Recently, 



further progress has been made to obtain the analytic finite temperature thermodynamics 



and quantum criticality of ID attractive fermions with strongly attractive interactions 
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29|. 



(l+l)-dimensional critical systems not only have global scale invariance but exhibit local 
scale invariance (conformal invariance) too. The conformal group in (l+l)-dimensions is 
infinite dimensional and completely determines the critical exponents and bulk correlation 



functions at criticality for gapless excitations 



30l |. Close to criticality, the dispersion relations 



for ID quantum systems are approximately linear. Conformal invariance predicts that the 
energy per unit length has a universal finite size scaling form E = Eq + A/L"^ where Eo is the 
ground state energy per unit length for the infinite system and A is a universal term. These 
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universality classes are characterized by the dimensionless number C 



contained in the term 



A) which is the central charge of the underlying Virasoro algebra |3ll. |32|. Affleck 3l| also 
showed that conformal invariance gives a universal form for the finite temperature effects 
on the free energy by replacing 1/L with T in the conformal map z = exp(27rtj/L). At 
the same time, Cardy 33|] showed that the two-point correlation function between primary 
fields can be directly derived from conformal mapping using transfer matrix techniques and 
expressed the conformal dimensions in terms of finite-size corrections to the energy spectrum. 
When C < 1, it takes on discrete values only i.e., C is cjuantized and hence the conformal 



dimensions are restricted to certain rational numbers [34\. On the other hand when C > 1, 
the critical exponents may depend continuously on the parameters of the model 35 1. 

The critical exponents for BA integrable models can be calculated via the quantum inverse 
scattering method (QISM) in terms of a function Z{X) called the dressed charge. In this way 
Bogoliubov et al. 36| obtained explicit expressions for the correlation functions for the Rose 
gas and the XXX and XXZ chains. Izergin et al. js^l considered the finite-size corrections 
to multicomponent BA systems and presented a formula for the dressed charge matrix Za^ 
which determines the critical exponents. They also showed that the integral equations for 
the dressed charge matrix depend only on the quantum i?-matrix of the model, which means 
that the universal class of critical exponents are described by the i?-matrix and the structure 
of the ground states of the integrable models. This universality pro per ty is a consequence of 
conformal invariance. Other models like the impenetrable Bose gas 



t — J model 



39| and the Hubbard model [40 
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the supersymmetric 



45| have also been considered in the 



context of the QISM approach. The Luttinger liquid is an alternative approach based on the 
fact that these models are certain realizations of the Gaussian model H] • Progress has 
also been made using Fredholm determinant representations of time-dependent temperature 
correlation functions for bosons and fermions in ID when c = oo 
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In this paper, we focus on the universal nature of ID repulsive spin-1/2 fermions in the 
frame work of the TBA formalism, including Luttinger physics and critical behaviour of 
correlation functions. In order to elucidate the significant features of spin-charge separation 
and critical exponents at quantum criticality, it is essential to analytically calculate the 
dressed energy potentials which encode the quantum and thermal fiuctuations of the spin 
and charge degrees of freedom in the critical regime. We investigate the low temperature 
thermodynamics of strongly repulsive spin-1/2 fermions in a small magnetic field and a 
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magnetic field close to the saturation field via the TBA method. We take an approximation 
to the TBA equations in the strong coupling regime, where the interacting strength c ^ 1. 
Thus the TBA equations are transformed into a new set of equations which can be solved 
using the Wiener-Hopf method. A comparison of the pressure and the entropy is made be- 
tween the application of two different integral expansions. These are Sommerfeld's lemma, 
which is valid for very low temperatures, and the polylogarithm function which is valid for 



finite temperatures. T 
theory prediction jsil 



le result from Sommerfeld expansion agrees with the conformal field 



32| . It is shown that the low energy physics can be described by a 



spin-charge separated theory of an effective Tomonaga-Luttinger liquid and antiferromag- 
netic SU{2) Heisenberg spin chain. A universal crossover from a relativistic dispersion to 
a nonrelativistic dispersion is determined by the exact thermodynamics extracted from the 
polylogarithm function. We also derive the explicit dressed charge matrix elements for the 
model in a weak external field H <^ 1 using the Wiener-Hopf method again, and also for 
the case that is close to the ferromagnetic state H — )■ He-, where He is the critical magnetic 
field. Various two-point and multi-point correlation functions at zero temperature are de- 
rived based on the expressions obtained from the dressed charge matrix. The leading terms 
and their critical exponents are given explicitly. Our results show that there is indeed no 
long-range order in this system. 

The spin-1/2 fermion model under consideration is the continuum limit of the ID Hubbard 
model (see, e.g., pp 45-49 of Ref. 6|), which has been widely studied jsl. In particular, the 
various correlation functions and scaling dimensions obtained here for the spin-1/2 fermion 
model have been derived for the ID Hubbard model, for interacting fermions and for a 
mixture of bosons and fermions using the dressed charge formalism [6|, 
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43| . Accordingly 



our results in the infinite coupling limit reduce to those obtained for the ID Hubbard model 
with an infinitely strong repulsion. Caution should be paid to the order of the limits T — )■ 
and infinite strong coupling [50]. Taking the T — )■ limit first, the correlations (for example, 
the one-particle correlation) show the scaling behavior of conformal field theory in the infinite 
strong coupling limit. However, taking the infinite strong coupling limit first, the correlations 
decay exponentially in the T — )■ limit. The two limits do not commute. Moreover, in the 
grand canonical ensemble, the Tomonaga-Luttinger liquid exists only in a certain region 
where the chemical potential is greater than the critical value. Below the critical chemical 
potential, the low temperature thermodynamics is that of an ideal gas in another regime. 
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see [51|. 



The corrections we have obtained in terms of strong but finite couphng are of impor- 
tance because of the experimental developments which enable access to the finitely strong 
coupling regime [szl . Our analytical 1/c order corrections to the critical exponents indicate 
an important signature - the critical exponents depend on the model parameters with cen- 
tral charge C > 1. In addition, our thermodynamical properties are valid for temperatures 
from T = to T ^ c^. Results of this kind are necessary to test critical phenomena and 
spin-charge separation theory in experiments with trapped fermionic atoms. Indeed, along 
with extending the known results for the thermodynamics and correlations, this is our main 
motivation here. 

This paper is set out as follows. In Section [IT] we introduce the model and present the 
TBA equations. The low temperature thermodynamics is derived in Section HTTl by expressing 
the TBA equations in the form of Wiener-Hopf integral equations. We solve the dressed 
charge equations for the model in the small field limit {H ^ 1) and in the limit where the 
field approaches the critical value {H — )■ He) in Section [IVl In Section |Vl we calculate the 
correlation functions of various operators in both limits. We then give a summary of our 
main results and concluding remarks in Section IVII Some detailed working and results are 
given in the Appendices. In Appendix |Al we derive the ground state thermodynamics and 
the critical field for the model. The Wiener-Hopf method is discussed in Appendix |Bl A 
more detailed derivation of the low temperature thermodynamics is given in Appendix O 
The leading terms of the zero temperature correlation functions are given in Appendices [Pl 
and El 

II. THE TBA EQUATIONS 

We consider a system of ID spin-1/2 fermions with delta-function interaction, with hamil- 
tonian 



Here = A^^ + A^^ is the total number of spin-up A^^- and spin-down fermions and 
M = (A^ — Ni)/2 is the magnetization, /i and H are the chemical potential and the 
magnetic field. In this paper, we exclusively consider the case of repulsive interaction for 





(1) 
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which c > 0. 

The ground state properties can be obtained using the BA solution 
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5J1 (see Ap- 



pendix [X] for a brief review). On the other hand, at finite temperatures T > 0, the physical 



quantities are described by the fo^ 
ferred to as the TBA equations 



lowing set of non-linear integral equations, which are re- 



59| . In the thermodynamic limit, their explicit form 



IS 



H 



e{k) = k'^ - ji T ^ a„ * In (l + e" 



4'n{k)/T 



) 



n=l 



0„(A) = nH- Ta^ * In (l + g-^^^)/^) + T T,^ * In (l + e-^-^")/^) 

m=l 

or equivalently 

e{k) = k^-n-TK* ln(l + 6"^^/^) - Ts * In (l + e^^^'^)/^) 
01 (A) = Ts * In (1 + e<^^(")/^) - Ts * In (l + 6"^^/^) 
0„(A) = Ts * In (1 + e<^"-i(^)/^) + Ts * In (l + e-^'^+i^^)/^) 

where the functions 0„(A) must satisfy the condition 

n— s>oo fl 

The functions a = a„(A), s = s(A), K = K{X) and Tnm = ^nm(A) are defined by 



(2) 

(3) 



(4) 
(5) 
(6) 



(7) 



a„(A) 



nc/2 



TT {ncy/A + A" 



1 



2ccosh(7rA/c) 



1 



27r / ^ 1 + e^^l^ 



a|„_„l(A) + 2a|„_^l+2(A) + . . . + 2a„+m-2(A) + a„+m(A) for n 7^ m; 

Tnm{X) = \ (8) 

2a2(A) + 2a4(A) + . . . 2a2„_2(A) + a.2n(A) for n = m. 

The asterisk * denotes the convolution / * g[x) = f{x — x')g{x')dx' . 

The bulk quantities are characterized by the solution to the TBA equations. For instance, 
the free energy per unit length F and the pressure P are given by 



rji /»00 

F = fm,-P, P=— ln(l + e-'^^'^''^)dk, 
where denotes the particle density. 



(9) 
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III. LOW-TEMPERATURE THERMODYNAMICS 

The most complicated part of the TBA equations is the string part (which characterizes 
the spin excitations) consisting of an infinite number of string functions 0„(A) (([3]) or ([S]) 
and ([n]))- These coupled nonlinear integral equations have not been solved in the most 
generic manner, but the obstacles can be overcome if we make certain assumptions for the 
parameters involved. Among them, one of the most crucial cases which we consider below, is 
the low-temperature limit T <^ 1. In this limit, the TBA equations reduce to a set of linearly 
coupled equations, which are easier to deal with. Moreover for the strong coupling regime 
c 1 in a weak magnetic field H <^ 1, we can solve the hnear integral equations analytically. 
Below we derive the analytical solutions in this physical regime: strong coupling c ^ 1, weak 
magnetic field H <^ 1 and low-temperature T ^ 1. The low-temperature thermodynamics 
for the generic case of c > and H < {He is the critical field) is derived in Appendix O 

We first observe from equation ([6]) that 0„,(A) > for n > 1, because s(A) > and 
ln(l -|- e'^"'^'^)/"^) > for A G M and n > I. This positivity condition implies that the function 
Tln(l + e-'^"^^)/^) for T ^ and n > 1. Therefore, in the low temperature limit 
T ^ 1, all the higher spin string functions drop off leaving only the function 0i(A) in the 
first set of TBA equations. The revised form is 



We now have to solve two coupled integral equations with only two unknown functions 
e{k) and 0i(A). Let us analyze them in the strong coupling regime c ^ 1, where the above 
equations are further simplified. Analyzing the dispersion e{k), we are able to rewrite the 



where the major contribution to the integral comes from a finite range {—ko, ko). As shown 
in Appendix lAl one notices that the points ±kQ correspond to the Fermi points in the charge 
Fermi sea. Thus equation ( ITTj) simplifies to 



6{k) = e~fi- — -Tai* ln(l + e--^^^*^)/^) 
0i(A) = H-Tai* ln(l + e'^^^)/^) + Ta^ * ln(l + e'^^^^)/^). 



(11) 



(10) 



term 55 1 




(12) 



0i(A) = H- 27rPai(A) + Tag * In (l + e-*'^^'>^^) . 



(13) 
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Using the Fourier transform, this equation and Eq. (fTOl) can also be expressed as 

£{k) = e -fi- 27TPK{k) -Ts* ln(l + e'^^^'^)/^) 

0i(A) = ^- 27rPs(A) + TK * ln(l + e^^^^^/^). (14) 

To proceed further, let us separate the function (f)i{X) into two parts: 

</.i(A) = </)f)(A) + rf)(A). (15) 

The first part ^['^■'(A) corresponds to the leading order term when T = 0, while the second 
part (j)^^\X) is the first order correction to the limit T — )■ 0. Analyzing the leading term 
H/2 — 27tPs{X) in Eq. (fT4|) . we find that (j)f'\X) should satisfy the linear integral equation 

0f HA) = y - 27rPos{X) + K * 0f )+(A) (16) 

where Pq denotes the pressure at T = 0. Here we have divided (j)f\X) into its positive and 
negative parts: 

^(o),u = <^(o)+(A) + 0(o)-(A), 

\(l>f\X) for |A| < Ao ^^^^ 
for |A| > Ao 

Note that the function (t)f\X) is nothing but the dressed energy es{X) ( 1A4I) denoting the 
energy of a spinon excitation with rapidity A, and the points ±Ao are the Fermi points (see 
Appendix \K\ for details). On the other hand, the function e{k) at T = corresponds to the 
dressed energy edk) flA4p describing the energy of charge excitation with momentum k. 

Substituting Eq. ( IT5|) into Eq. (fT4|) . and subtracting Eq. f|T6|) from the resulting equation 
gives 



0i')(A) = -27r(P - P,)s{X) + [ ir(A - /i) (Tin f 1 + e(^i°'('')+*^''(^))/^) - 0f ^+(/i)) rf/i 
+ / Tir(A-/i)ln(l + e<'('^)+'^i"('^»/^)rf/i. (18) 

An iteration procedure shows that 0[^^(A) = o(T). Thus one sees that 

0f ^A) ^ -27r(P - Po)s(A) + Ek[X) + [ K{X- fi)<p'l\fi)dfi (19) 



where we denote 

Ek{X) = TK*\n{l + e^<^^°' W'/^) . (20) 

In the hmit T — )■ 0, the major contributions towards Ek{X) come from the regions near 
±Ao. Hence we expand 0[°'*(A) around A = ±Ao 

</)f)(A) = t(A-Ao)+0((A-Ao)^) (21) 

where t = d(j)f'\X) / dX . Then we find 

A=Ao 

Ek{X) I K{X-^i) ln(l + e-*l'^-^°l/^)ci/i 

J \fi~\o\<e 

r)T-i2 poo 

~ [K{\ - Ao) + K{X + Ao)] / In (1 + e"") du 
t Jo 

= fIl[K{X-Xo) + KiX + X,)]. (22) 
or 

Therefore we obtain a hnear integral equation which determines (j)^^\X), namely 

4'\X) ^ -27r(P-Po)s(A) + ^[K(A-Ao)+ir(A + Ao)]+ / KiX- fi)(t>? Wfi. (23) 

or Jm>Xo 

In completely the same way, one finds that the low temperature behavior of e{k) (IT^ is 

described by 

e{k) ^ k^-fx-2nPKik)-^^[s{k-Xo) + sik + Xo)]- [ s{k~X) Uf)(A) + rf^(A) 

6^ ^|A|>An L 



dX. 

(24) 

a. firoil and Eq. f l2^ can be solved for H 1 via the Wiener-Hopf technique as in 
. |20| ■ For convenience, let us introduce the functions 

l/«(A):=0S'=)(A + Ao) (A: = 0,1). (25) 

By definition (see Eq. f|T7|) ). the Fermi points ±Ao are determined by the condition 

0f)(Ao) = y(°)(O)=O. (26) 

Applying the iterative procedure to Eq. (fT6|) . one sees t/*-°^(0) = = H/2 — 2txPqs{Xq) + 
0{H). Solving this equation, one finds that Aq ~ — Inif for H <^1. Because K{X) rapidly 
decreases with A > 0, we may solve the integral equations (fT6!) and (!23|) by expanding 



^'\X) = Y.yn\^) (^ = 0,1). (27) 



y 

n=0 
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yn\X) obeys the integral equation 



(A) = gl^KX) + I K{X- 



(28) 



where the driving terms gn\X) and gn'{X) in the hmit of T <^ 1 are exphcitly given by 



gl^\X) = ^j--2nPos{X + Xo) 
gl,'\X) = -2niP-Po)s{X + Xo) + 



H 



6t 



■KiX) 



POO 

(A) = / KiX + fi + 2Ao)yJi (/i)t//i 

^0 



(n > 1). 



(29) 



The above equations are the so-called Wiener-Hopf type integral equations. In Ap- 
pendix [HI a method to solve the Wiener-Hopf type integral equations is given. Let us 
decompose y^^\(jj) into a sum of two parts, i.e., if^\uj) = y^\uj) +yL\uj), where ]/f!\uj) 
(^■'(w)) is analytic in the upper (lower) half plane (see Eq. (1B3P in Appendix [B]). For the 
leading terms yg^''(a;), we obtain (see Eqs flB20p and flB23p ) 

27riPoG'_(-7ri/c)e-'^^«/'=" 



-le] 



2{u + ie) 
{G+iu) - 1) 



+ 0{H' 

c uj + iri/c 

2m{P - Po) G+(u;)G_(-7ri/c)e-"^«/'= 



(30) 



Qt ^ ' ' c UJ + 111/0 

where G±{uj) is defined by Eq. (1B12|) . Using the formula as in (IBlOp . and combining the 
above equation with the condition (126|) . we determine the leading term of the Fermi points 
Ao to be 



An--ln(^ 



TT 



H 



Hn 



4:7iPoG-{-7ii/c) AnPo 



'TT' 



- = J-H, + 0[-\ (31) 



1 



cG'_(0) c V 2e V 2e' 

where we have used the formulae (IB5I) and flB13l) . To derive the last equality in the second 
equation we used the property (see Appendix [A]) 



(32) 



where is the critical magnetic field, where all fermion spins point up (see 57| or Eq. (IA9 
in Appendix El) . with He ~ 8nj?7r^/3c. 
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From the relation ( 125|) . the integral in can be evaluated as follows: 

9„— -Ao r(X) 

s{k - \)<pf\X)d\ ^ \ e-^^\f{\)d\ 

|A|>Ao C Jo 



C ./n .Loo 



ze-^^« r dco ^k) 



, (33) 

where the relation s{k ± |A + Ao|) ~ g-'r(A+Ao)/c^^ used in the first line. Substitution of 
Eq. ([33]) together with Eqs. 031]), (130]) and flBTSj) into Eq. ^ yields 

4t)^t^-,-2.PK(k)-g^-^^. (34) 

Inserting the relations 

Kik) ^ — , y',{0) = - hm u'yS'^iu) = ^ (35) 

and using P Pqi "we arrive at 

,2 2Pln2 cH^ cT^ ^ 

_ — . — ^e-A (36) 

where A := yU+2Pln2/c+cif^/47r^P+cT^/6P. We would like to address that the calculation 
of the dressed energy potential ( 15B]) is essential for catching spin-charge separation signature 
and quantum criticality at low temperatures. The terms in the function A show an important 
implementation of spin density and charge density fluctuations that reveals a physical origin 
of spin-charge separation. This is clearly seen from the following low temperature and finite 
temperature thermodynamics. In view of this validity of catching this universal nature, we 
see that the result (!36l) is also helpful to numerics. Using this function A, we then perform 
integration by parts on the pressure of the system (j9]) to get 

where ljis{z) is polylogarithm function defined by 

k=l 



12 



For fixed particle density, the chemical potential is determined by solving the equation 
rtc = dP/ dfi. By inserting Eq. fl37|) into Eq. (Q , the low-temperature thermodynamics in the 
region ^ 1 and c ^ 1 is completely determined. To capture universal features of the low- 
temperature thermodynamics, we further expand Eq. (13 7p by making use of Sommerfeld's 



lemma (see 



55 



56( 1 for instance): 



Furthermore, we use the relation Uc = dP/dfi and repeatedly iterate the terms to elimi- 
nate those with orders higher than T^, and 1/c. After some lengthy algebra, we finally 
obtain the chemical potential 

2 2/^. 161n2\ ^iH^ 31n2\ ^T^ / 31n2\ 

where 7 denotes 7 := c/ric- Substituting (HOj) into the expression for the pressure (!39ll and 
then iterating the pressure itself, we obtain 



2 „ 61n2\ 97/^2 / 41n2\ 37r2 / 41n2\ ^ ^ 



3 ^ V 7 y Stt^tIc \y ' 7 y Att'^Uc \ 1 J ' Sric 



The free energy of the system defined by Eq. flC14p is given by 

I 2 sf. 4:\n2\ , Q\xi2\ 7T2 Q\n2\ 



F ^ -vrX 1 -r— 1 + - 1 + • (42) 



From this free energy ( 142|) . we have the susceptibility 

<92F 37 / 61n2\ 

This susceptibility can be possibly tested in an trapped ID Fermi gas of cold atom. We 
rewrite the pressure ( 1371) 

m „3, 



where 



^ /i ln2 /T / 1 H 

X = —\ — Li3 — ' 



T v^Veo ; 27rfT2 /ZLi3(-e^) 

,„Li| (-e^) 



(45) 
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In dimensionless units, i.e. h = H/eq and ^ = h/sq with Eq = |^c^, the susceptibihty is 



X c 37 Ai + 61n2, (45) 



given by 

which agrees with the field theory prediction x^s = 0/11. Here 9 = 1/2. The holon velocity 
Vc and spinon velocity Vg at H = are given in (1C16I) in Appendix O In the above equation, 
the density is given by 



, n II 
n : = ^ — 



3 fhyf^^3^\n2\^8^\n2 



and the dimensionless interaction strength is given by 

1 3 fh 



7 TT Svr^ \/i 



(47) 



(48) 



In Fig. [H we plot the susceptibility for different values of the chemical potentials at low 
temperatures. In contrast, in Chapter 13 of Ref. pj, the thermal and magnetic properties 
of the ID Hubbard model are plotted by using the result obtained for the thermodynamics 
of the ID Hubbard model by numerically solving nonlinear integral equations (NLIE). The 
thermodynamics for the Hubbard model can be calculated by the quantum transfer matrix 
method for all temperatures. Here we have obtained an explicit low temperature expansion 
for the fermion model. The result (jSj) for the pressure contains the spin density and charge 
thermal potentials at criticality and may thus possibly be used to test the spin and charge 
velocities in experiments with ultracold atomic fermions in a ID harmonic trap. 

By substitution of the holon velocity Vc and spinon velocity Vs at H = (see (1C16P in 
Appendix[C]), the free energy P2|) suggests the universal low temperature form of spin-charge 
separation theory, namely 

(49) 



6 

where Eq denotes the ground state energy density and C = 1. The above behavior can also 
be derived for generic c and H (see Appendix [C]) . This expression corresponds to two central 
charge C = 1 conformal field theories 3l|]- This universal nature of ( H9|) means that the 
low-lying excitations are decoupled into two massless degrees of freedom which are described 
by two Gaussian theories. However, if the excitations involve highly excited states, these 
theories break down. In Fig. [2] we compare the exact thermodynamics with the predictions 
of conformal field theory and show the breakdown of the two Gaussian theories at higher 
temperatures. 
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FIG. 1: Susceptibility vs temperature for different values of chemical potentials /i = 0.01, 0.03, 0.05 
and h = 0.001. At T = 0, the susceptibility values for different values of chemical potentials are 
consistent with the field theory prediction x^s = O/tt- This possibly gives a way to test effective 
spin velocities of the spin-1/2 ultracold atoms with a repulsive delta- function interaction. 

IV. DRESSED CHARGE 

The excitation spectrum for repulsive spin-1/2 fermions is gapless for any magnetic field 
strength H. As a consequence, the asymptotic behavior of the correlation functions of the 
model can be described by conformal field theory (CFT) 3ll-l33|. CFT relates the critical 
exponents of the correlation functions to the finite-size corrections in the energy spectrum. 
The basic tool used to determine the critical exponents from the finite-size corrections is the 
dressed charge matrix. For this model, it is a 2 x 2 matrix which couples the charge and 
spin degrees of freedom together and at the same time governs the excitations of charge and 
spin waves near the Fermi surface. 

The dressed charge matrix of this system is explicitly given by 

Z = I ^''^'"'^ I (50) 
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Entropy Polylog 

Entropy Sommerfeld 

Pressure Polylog 

Pressure Sommerfeld 



_l I I L_ 



_l I I L 



J I I L_ 



_] I I L_ 



_] I I L_ 



0.4 0.6 
Temperature T 



FIG. 2: (Color online) Pressure and entropy vs temperature T. Solid lines show the pressure (j4ip 
and entropy obtained from the free energy (j49p : dashed lines show pressure and entropy obtained 
from the polylogarithm function (j37p which gives the precise thermodynamics for temperatures 



2 fc2 2 

below the Fermi temperature ksTp = Ep = ^ in the strong coupling limit. At low temper- 
atures, they are indeed consistent with the field theory predictions. The deviation of the entropy 
from the linear temperature dependence marks a breakdown of the two massless field theories. 



while the integral equations of its elements are given by 



Zccik) = 1+ r a,{k-\)Z,s{\)d\ 

J-Xo 



ai{X- k)Zccik)dk - / a2{\ - lj)Zcs{lj)diJi 

-ko J—Xo 

/Xo 
ai{k - \)Zss{\)d\ 
Xo 

/ko pXo 
ai{X - k)Zscik)dk - / a2(A - 
'ko J -Xo 

This set of equations are in turn made up of two coupled sets of equations. Equations (|5T|) 
and (152!) can be treated separately from equations (!53|) and (15^ . The following relations 



(51) 
(52) 
(53) 
(54) 
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FIG. 3: (Color online) This figure shows the dressed charges Zcdko), Zcs(Ao), Zsdko) and Zss(Ao) 
as a function of the external field H/Hc- The dressed charges are plotted for different values of 
Uc/c which is the inverse of the interaction parameter 7. These curves are plotted by numerically 
solving Eqs. ([5T]) -(l5T 



are useful for further calculation 



-ko 27r 
2tx 



dk 



dk 



Zcs{\) 

27r 



d\ 



Ao 



ps{\)d\ = ni 



Ao 



Pc{k)dk = ric 



(55) 
(56) 



where ni is the density of down-spin fermions. To derive these relations, we first multiply 
each set of density equations with the dressed charge equations and integrate them. While 
making use of the fact that the kernels are symmetric, we then subtract one equation from 
the other to eliminate the same terms. FIG. |3] shows a depiction of the numerical solutions 
to Eqs. dSID-dSr 
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A. The limit H <^ 1 



When H = 0, the ground state of the system is antiferromagnetic. We now solve the 
dressed charge equations under a small magnetic field by the Wiener-Hopf technique. Con- 
sidering only terms of up to order 1/c in the strong coupling limit, equation can be 
written as 



//Co /-AO 
Zsc{k)dk - I a2(A - n)ZssifJ^)dn 
'ko J -^0 

Ao 

where the property fl55|) was used. Applying the Fourier transform, we obtain 
ZssW = \ + 2vm4s(A) + I K{\ - ^)Zss{fi)dfi. 



(57) 



(58) 



This equation is the same as Eq. ( |T6l) other than the driving term, and hence a similar pro- 
cedure as introduced in Sectionlllllis applicable. Introducing the function y{X) = Zss{X + Xq) 
and expanding it as y{X) = Yl'^=o VnW, "we obtain the integral equation (see the correspond- 
ing equation fl28|) in Sec. IIIII) 

/■oo 

(59) 



VnW = fl'n(A) + / K{X - i2)yn{ij)dii 
Jo 

where gn{X) is given by 

9oW = ^ + 27™|s(A + Ao), 

POO 

g„{X.) = I f('(A + /i + 2Ao)y„-i(/i)(i/i {n>l 
Jo 

For n = 0, setting a = 1/2 and b = 27m in Eq. ( ]B20p . one has 

'iG_(-ie) 27rm^G'_(-7ri/c)e-^^»/'= 



+ 0{H' 



(60) 



(61) 



2{uj + it) c{uj + iri/c) 

Combining this result with Eq. (IBlOp . and finally substituting the relations (IBlSp . (1B5P and 
( l3T|) . one obtains the first order contribution to Zss{Xo): 



*(0) ^ (l 



H 



0{H^). 



(62) 



To obtain the second order correction to y(0) = Zss{Xq), we must consider the contribu- 
tion of i/i(0). The Fourier transform of 5'i(A) is given by 

1 



-uj)K{uj) 



1 - 



G+(w)G_(w) 



(63) 
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Here we have used the decomposition of the kernel (IBSp . As demonstrated in Appendix. [Bl 
let us decompose ^i{uj)G-{u) into the two parts $±(a;) which are analytic in the upper and 
lower half planes, respectively: 'gi{u)G^{u) = $+(w) + ^^{u). Now $+(0;) is given by 



-dx 



(64) 



2ni X — u — ie G+{x) 

where e is a small positive constant. Note that the function has a branch cut along 

the negative imaginary axis. Deforming the integration contour to avoid the branch cut, we 
have 

^' " 1 1 



1 

27ri jQ 
1 

(27r)2? Jo 
2 



X — lUJ 



X 



lUJ 



{2n 
1 

1 



~ e-^^^^'yo+iix) 



3 

2 Jo 

00 g-2Aox 



C 



a; — 20; 
c 

+ 



V27re 



— e) G^ 

CX / -1 

27r _ / 1 

2 



-ix + e) 



dx 



/ CX \ ^ ^ / i CX 



27r 



e 2 



— tca^ \ 

e 2 jctx 



27re. 
+ 0(x) 1 dx 



1 CX \ / CX ^ 

r ( I sin ( — \ dx 

2 27r / V 2 



(65) 



-iw V4A/27rAo 

Note that, from the third to the fourth step in the above equation, we have used the 
fact that the integrand rapidly decreases for x > because Aq ^ 1, and hence the integral 
can be approximated by expanding the terms other than exp(— 2Aox) around x = 0. By a 
relation similar to Eq. flB9p . yi+(u;) is expressed as 

GAu\ 



-loj V4v27rAo 
Insertion of the relation flBlOp and Eq. flB5l) yields 



«4 



(66) 



Therefore from 



(1^ and (IHTj) . we finally obtain the expression 
4ni f H \ 1 



(67) 



1 
71 



1 + 



+ 



+ 



{\nHo/Hy 



(68) 



c \HJ ' 4\n{Ho/H)_ 

Next we evaluate the dressed charge Zsc{ko). The Fourier transforms of (153|1 and ( 158|) 
give 



Z,,{k) = i + 27™;K(A;) - / 



s(A;-A)Z,,(A)rfA. 



(69) 



\li\>Xo 



19 



Applying the same procedure as used in the derivation of Eq. (l33i) . we immediately obtain 

1 _ . 2e-^^V 



+ 



2ni In 2 2 



c 



HMHo/H) 



2 c 

where we have used the relation fl35l) . 

Repeating this whole process to evaluate Zcs{\o) and Zcc(/i;o), "we find that 



(70) 



H 



and 



^cc(A;o) 



1 + 



272^, In 2 4nf. 



H,i\n{H,/H)f 



+ 



Hl{\n{H,/H)f 



(71) 



(72) 



The down-spin density can be explicitly written in terms of the external magnetic field 
H by evaluating Eq. ( 155|) . Using the property that Zgcik) ^ Zsciko) + 0(l/c^) for c ^ 1 
and fco ^ uric/ {I + 2 In 2/7) (see Eq. flC12l) ). we find that 



4 /if 



(73) 



By substituting the expression ( 1731) into Eqs. ( l68i) and ( ITOi) . the dressed charges in the strong 
coupling regime c ^ 1 and a weak magnetic field if ^ 1 are explicitly determined in terms 
of the fixed particle density and the external magnetic field H. 



B. The limit ^ for 7 » l/^l- H/Hc 

When H > He, the ground state of the system is ferromagnetic. Correspondingly the 
Fermi point Aq becomes zero. Before solving the dressed charged matrix for H approaching 
the critical field He from below, we have to know how Aq behaves in this vicinity. The spin 
part of the TBA equation at T = (i.e., (j)f'\X) in Eq. f|T6|) or equivalently the dressed 
energy es{X) in Eq. ( ]A4I) ) is approximately 

ai{\ — k)ec{k)dk — I a2{X — fi)es{fJ,)dfi 

^H- 271 Poa, (A) - 2 Aoa2 ( A)£, (0) (74) 

which is derived by approximating the first integral as in (fT2|) and expanding the second 
integral around Aq ~ 0. Let us assume 7 ^ 1/ ^^1 — H/ He ^1. In this region, one sees 
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from Eq. (1A9P that Eq. (17il) further reduces to 



r(\^-n 2Ao 5,(0) 

"^^^^ - ^ " TT4AW " ^i + aW ^ ^ 

With this result, we can exphcitly find ^^(O) to be 

i + zAo/vrc 

The last step to derive an explicit expression for Aq is to use the condition ^^(iAo) = 0. 
This gives 

TT 2Ao Hc-H ^ 

1 + 4A2/c2 ^ vrc (1 + A2/c2)(l + 2Ao/7rc) ' ^ ^ 

Multiplying the equation with each denominator and then ignoring the terms of order 
0(Aq/c^) or higher yields the equation 

H-H^+^-^H-^H^ = 0. (78) 

After rearranging the terms and using the fact that H He, we arrive at the result 



^««yi-f (79) 



4l| and the XXZ 



which is also similar to the result obtained for the ID Hubbard model 
Heisenberg chain jsGj. 

With this expression, we can evaluate the dressed charge matrix explicitly in terms of H. 
For Z,s{Xo), from ([52]) 

Z,,(A) ^ 27ia,i\)n, - 2Aoa2(A)Z,,(0). (80) 



Here an approximation similar to the above and the property (!56ll have been used. Solving 
for Zcs{0) and neglecting terms of order 0(l/c^), we have Zcs{0) ^ i/'j. Substituting this 
together with (I79|) imphes that 



2o.(A„)«i(l-yi-|j. (81) 

Zcc{ko) is easily obtained by substituting the above results into the approximation form 
Zccih) ^ 1 + 2Aoai(A;o)^cs(0) of Eq. The result reads 



Z,,{ko)^l + —Jl-^. (82) 
7r7 V He 
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The remaining dressed charges Zss{Xo) and Zsc{ko) are also evaluated by similar calcula- 
tions. For Zss(Ao), Eq. fl5^ becomes ZssW = 1 + 2'Kniai{\) — 2XQa2{X)Zss{0), and hence 
Zssi^o) = 1 + 4n|/c — 2Ao/(vrc). The density of down-spin fermions is evaluated by 
applying the same method as in the above to Eq. f lAip . with result 



= / p,{X)dX ^ 2AoP.(0) ^ -^Jl - — . (83) 



This in turn gives 



Likewise, 



I H 8 H 

Z„(A„)«l--^l-- + -^l--. (84) 



Z.o(*o)«^\/l-|j. (85) 

In FIG. m we compare the numerical solutions to the leading order solutions for the 
dressed charges in both limits H <^ 1 and H — t- He- It shows good agreement between 
our leading order solutions and the numerical solutions at points not far from H = and 
H = H,. 

V. CORRELATION FUNCTIONS 

We turn now to the calculation of the long distance asymptotics of various correlation 
functions and scaling dimensions which have been obtained for the ID Hubbard model in 
te™. of the dressed charge formal., fl U Q. Our results exte.d these results into tire 



strong but finite coupling regime for the spin- 1/2 repulsive fermion model. 

The spin-1/2 repulsive fermion model is gapless and thus critical at zero temperature. At 
T = 0, the correlation functions decay as some power of distance governed by the critical 
exponent which we shall denote conventionally by 6. For T > the decay is exponential. It 
was shown that conformal invariance leads to universality classes of critical theories that are 



related to the central charge C related to the underlying Virasoro algebra [30l]. The critical 
behavior of the model under consideration is described by the direct product of two Virasoro 
algebras, one characterizing the spin degree and the other characterizing the charge degree. 
Both Virasoro algebras have central charge C = 1. 

The general two-point correlation function for primary fields ip with conformal dimensions 



22 




N 



0.4 0.6 

H/Hc 




0.4 0.6 

H/Hc 




0.4 0.6 

H/Hc 




0.4 0.6 

H/Hc 



FIG. 4: (Color online) This figure shows a comparison between the numerical solutions (solid lines) 
and the leading order corrections to the dressed charges Zcc{ko), Zcs(Ao), Zsc{ko) and Zssi^o) in 
the hmits H <^ 1 and H ^ He for ric/c = 0.001. 



at T = and T > are given by 



and 



{x — ivcty^'= (x + ivcty^'' {x — iVsty^'' (x + ivst) 



2A: 



(86) 



X 



X 



exp{-2iDc{kp^ + kFi)x) exp{-2iDskFix) 



t>c sinh(7rT(x — iVct)/Vc) J \t>c sinh(7rT(x + iVct)/Vc 

2A^ 



vrT 



TT 



T 



TT 



T 



2A: 



2A.: 



(87) 



ysSmh.{TTT{x — ivst)/vs) J \Vs smh.{TrT {x + iVst)/vs 
where ^F^.t ^ire the Fermi momenta, < x < L and — oo < t < oo is Euchdean time. The 
conformal dimensions of the fields can be written in terms of the elements of the dressed 
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charge matrix as 

2Af{AN,,,, Ar± , D,,) = + Z^^D, ± ^^^^^^^f^) ' + (88) 

and 

2A±(A7V,_„ 7V± , i),^,) = (^Ze.L'c + Z.^D, ± ^ j + 27V±. (89) 

The non-negative integers AN a, and the parameter where a = c,s represent the 
three types of low-lying excitations. Here AN^ denotes the change in the number of down- 
spin fermions. characterizes particle-hole excitations where A^^ i^a) number of 
occupancies that a particle at the right (left) Fermi level jumps to. also enumerates 
the descendent fields for the primary fields (p. And lastly, Da represents fermions that are 
backscattered from one Fermi point to the other. They are restricted by the condition 

AN, + AN, , , AN^ , , 

Dc= (modi), Ds=^ (modi). (90) 

We want to find the asymptotic behavior of the general two-point correlation functions 
for the operators 0{x, t), namely {0{x, t)O^(0, 0)). The operators can be written as a linear 
combination of primary fields with conformal dimensions A^^ and their descendent fields. 
Noting that the correlation functions for fields with different conformal dimensions are zero, 
we can express the correlation functions at T = and T > respectively as 

{0{x, t)Ot(0, 0)) = y A{AN, N^,, ,) exp(-2zD.(fc^, + k,,)x) eM-^^Dskn^) 

(91) 



and 



(O(,x,t)Ot(0,0))T 

= J2 M^Nc,s, N^s, Dc,s) eM-2iDc(kFt + kFx)x) exp(-2iL>,A;F;x) 

n 

ttT \ ( T^T ^ 



VcSinh(7rT(x — ivct)/vc) J \WcSinh(7rT(x -|- ivct)/vc) 



X 

X ( — ) I — ) (92) 

\Vs smh.{7rT (x — ivst)/vs) J \Vs smh.{7rT{x + ivst)/vs) J 

where n denotes the set of quantum numbers 

n = (A7V,_„7V±,L»,^,) (93) 
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"or the form 



33| 



which are determined by the condition given in ( 190|) and the selection rules 
factors while performing a spectral decomposition of the correlation functions 

Let us consider the correlation functions of operators which are written in terms of the 
field operators ilj„{x,t) where a =t,i- They obey the canonical commutation relations 

{tp^{x,t),ijl{x,t')} = 6„„>6{x-x') 

{Mx,t),Mx,t')} = {iljlix,t),tljlix,t')} = 0. (94) 

Here we consider the following correlation functions: 

(i) One particle Green's function: 

G^{x,t) = {Mx,t)^l{0,0)). (95) 

(ii) Charge density correlation function: 

Gnn{x,t) = {n{x,t)n{0,0)) (96) 

where 

n{x, t) = ni-(x, t) + ni{x, t), n„{x, t) = ipl.{x, t)ipa{x, t). (97) 

(iii) Longitudinal spin-spin correlation function: 

G'{x,t) = {S'{x,t)S'{0,0)) (98) 

where 

S%x,t) = ^{n^{x,t)-n^{x,t)). (99) 

(iv) Transverse spin-spin correlation function: 

G^{x, t) = (S+(x, t)S-{0, 0)) (100) 

where 

S+{x, t) = ^lix, t)ij^{x, t), S-{x, t) = ^lix, t)^^{x, t). (101) 

(v) Pair correlation function: 

G,{x,t) = (V'4(x,t)V't(a;,^)W(0'0)^!(0'0))- (102) 
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For each of the correlation functions considered above, the values of n are given by 





[x,t) 


(AiV,= 


l,ANs 


= 0,D^eZ+ 1/2, DseZ+ 1/2) 


Gi 

4- 


[x,t) 


{ANc = 


1, AiV, 


= i,Dce z,Ds e z + i/2) 




t) 


(AiV,= 


0,AiV, 


= 0, G Z, G Z) 


QZ 


[x,t) 


(AiV,= 


0,AiV, 


= 0, G Z, G Z) 




^x, t) 


(AiV,= 


0,AiV, 


= 1,^)^ G Z + l/2,Ds G Z) 




[x,t) 


(AiVe = 


2,ANs 


= 1,^)^ eZ + l/2,Ds G Z). 



with iV^±, G Z>o for every case. The explicit results for these correlation functions for H <^ 1 
and H ^ He are given in Appendices [D] and [El which include the order of I/7 corrections 
in the critical exponents. 



VI. CONCLUSION 



We have derived the low temperature thermodynamics and long distance asymptotics 
of correlation functions for the spin-1/2 repulsive delta-function interacting Fermi gas with 
an external field by means of the thermodynamic Bethe ansatz method and dressed charge 
formalism. With the help of Wiener- Hopf techniques we have calculated the low temperature 
free energy and thermodynamics and found that the low energy physics can be described by a 
spin-charge separated theory of a Tomonaga-Luttinger liquid and an antiferromagnetic spin 
Heisenberg chain. The dressed charge equations have been solved analytically for a small 
external field if — )■ and a large external field H ^ He using the Wiener-Hopf method. We 
have also calculated the conformal dimensions for many correlation functions including the 
one particle Green's function, the charge density correlation function and pairing correlation, 
as given in Section |Vl 

In particular, the explicit form of the critical exponents which we have obtained in terms 
of the external magnetic field and the interaction strength up to 1/c corrections extends 
the known results obtained for the ID Hubbard model in the infinite coupling limit [gI. l40|. 
They provide insight into understanding the critical behaviour of interacting fermions in ID. 
The result for the free energy at low temperature shows a universal signature of Tomonaga- 
Luttinger liquids where the leading low temperature contributions are solely dependent on 
the spin and charge velocities. It is to be expected that this universal nature can be tested 
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via the finite temperature density profiles of tlie repulsive Fermi gas in an harmonic trap. 
This opens a way to experimentally observe how the low temperature thermodynamics of a 
ID many-body system naturally separates into two free Gaussian field theories. 

We have also presented results for the low temperature thermodynamics which extend 
beyond the range covered by spin-charge separation theory. As effective as it is, the Wiener- 
Hopf method does not allow the full derivation of the equation of state for temperatures 
beyond the Tomonaga-Luttinger liquid regime. This is because the Tomonaga-Luttinger low 
temperature physics does not contain enough information on thermal fluctuations necessary 
to describe the quantum critical regime. This restricts access to quantum criticality in 
the whole physical regime. However, as we have demonstrated, the polylog formalism is 
suitable for the study of low temperature and strong coupling in the fi — H phase plane 
for weak magnetic fleld. The pressure we obtained captures the essential spin density and 
charge density fluctuations at criticality and may possibly be used to test the spin-charge 
separation theory in experiment with ultracold atomic fermions. This points to the further 
study of quantum criticality in ID interacting Fermi gases with repulsive interactioii, as has 
been done recently for attractive interaction and for the Lieb-Liniger Bose gas 28|, |29 |. 

This work has been supported by the Australian Research Council. The authors thank 
Professor Tin-Lun Ho and Professor Rafael I. Nepomechie for helpful discussions. 
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Appendix A: The ground state properties 



Here we briefly describe the ground state properties of spin-1/2 fermions witli repulsive 
interaction. The full spectrum of the Hamiltonian can be obtained by the BA method 

~i n 

531 . I54j |. In the thermodynamic limit, the ground state properties are characterized by 
two Fermi seas made up of charges and (down) spins. The distribution functions Pc{k) of 
charges with holon momentum k, and Ps{k) of down spins with spinon rapidity A are written 
as integral equations, 



1 



Ao 



Pcik) = — + ai{k - X)ps{X)dX 
2^ J-Xo 

/ko pXo 
ai{X - k)pc{k)dk - a2{X - p)ps{n)dn (Al) 

where the function a„(x) is given by ([8]) and ifeg and ±Ao correspond to the Fermi points. 
The density of the fermions Uc = + n-i {ricr denotes the density of spin-cr fermions) and 
the density of the down-spin fermions are respectively given by 

ko rXo 

Pc{k)dk = nc, / PsiX)dX = n^. (A2) 

■ko J ~Xo 

The ground state energy per unit length (denoted by Eq) is 

/ko rXo 
{k'^ - p- H/2)p^{k)dk + H psiX)dX (A3) 
'ko J -Xo 

where Pq is the pressure at zero-temperature (see ([9]) for finite temperature). The ground 
state properties are also described in terms of the charge dressed energy edk) and the spin 
dressed energy es{X) as 



e^{k) = e -p- H/2 + ! ' ai{k - X)es{X)dX 

J~Xo 

/ko rXo 
ai{X — k)ec{k)dk — / a2{X — p)es{p)dp. (A4) 
-ko J-Xo 

The above dressed energies define the energy bands. The ground state corresponds to 
the fillings of e^k) < and es{X) < 0. Thus the Fermi points ±ko and ±Ao are given by 
the conditions 

ec{±ko) = 0, £.(±Ao) = 0. (A5) 
Using the above dressed energies, Eq is written as 

^0 - pnc = -Po = Tr / £c{k)dk. (A6) 
2vr J-ko 
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One can immediately realize that the dressed energies edk) and ^^(A) respectively correspond 
to 6{k) and 0i(A) in the TBA equations ([2]) and ([3]) in the limit T — )■ 0. 

Let us evaluate the value of the critical field H = He- At this point, the density of 
down-spin electrons is zero {n^ = 0) i.e., Aq = 0. Therefore the expressions (lAip and f lA4p 
are significantly simplified. Inserting Pc{k) = l/(27r) into (IA2p . one has ko = irnc- The 
Fermi point ko is also calculated by edk) = k"^ — ji — Hc/2 (]A4p and the condition (lASp . i.e. 
kg = fi + Hc/2. Substituting these expressions into ^^(A) and using the condition (lASp . one 
finally arrives at 

H^= ( ^ + 27^2 ) tan"^ / ] _ ^j^j^ 



,27r 

The pressure Pq at H = is 



Po = '^n^nl (AJ 



3 

In the strong coupling limit 1, and Pq are given by 

^^^V^-^J' + ^^^^ 

Appendix B: The Wiener-Hopf method 

In this appendix, we briefiy review how to solve the integral equations appearing in the 
main text (see Eq. f l28|) for instance) by using the Wiener-Hopf method. 
Consider a Wiener-Hopf integral equation 

y{X) = g{X)+ / ir(A-^)y(/i)ci/x (Bl) 
Jo 

which determines an unknown function y{\) where K{X) is defined by Eq. (|8]), and the 
driving term g{X) defined on the entire real axis is assumed to be a known function. By 
Fourier transforming Eq. fIBip . we obtain 

y+M = (l - K{uj)y' {g{uj) - y_{uj)) (B2) 

where the functions ^±(0;) are defined by 

/oo 
e{±\)y{X)e'^-d\ (B3) 
-00 

with ^(A) representing the Heaviside step function. The function K{uj) can be explicitly 
written as 

^ -c\i^\/2 

K{u) = , , , , , (B4) 

^ ' 2cosh(c|a;|/2) ^ ' 
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Eqs. (1B3P denote a decomposition of y{ou) into the sum of functions analytic on the upper 
and lower half-planes, respectively, with y{uj) = y+{<jj) + y-{(jj). Hereafter we assume that 

y±(oo) = ?(oo) = 0. 

To solve the equation, we factorize the term 1 — Kiu)) as 

{l-K{oj)y^ = G+{oj)G_{oj), limG±(cj) = l (B5) 

where G'+(w) {G-{uj)) is a function which is analytic and nonzero in the upper (lower) 
half-plane. Since K{uj) = K{—u), one finds G^{u)/G^{—uj) = G^{—u)/G^{uj). Therefore 
G'+(co')/G_(— w) is a bounded entire function. Liouville's theorem and the asymptotics of 
G±{uj) yield 

G+{u) = G-{-u). (B6) 
Using the factorization equation (1B5I) . Eq. (1B2I) becomes 

- ^^co) = $_M - G-Mr (c^) (B7) 

where ^±{uj) are functions assumed to be analytic and bounded in the upper and lower half 
planes, respectively, 

G4uj)g{uj) = $+(0;) + (B8) 

FromEq. (IB7p . the function y+(tj)/G+(tj)—$+(a;) is a bounded entire function, and hence 
is a constant according to Liouville's theorem. Considering the asymptotics, we obtain 

y^iu) = G+H (<I>+H - $+(00)) , y4iu) = - <l>_(oo)) . (B9) 

A formula useful for practical calculations is 

yiO) = ^ r yi^y'^'dcu = ^ [ y4u)e-'-^dcu = -i lim ujy+{uj) (BIO) 

where G~ is a semi-circular path on the lower half-plane. We have used the fact that the 
sum of residues in the lower half-plane is equal to the residue at the point at infinity. 

Next we will determine the explicit forms of y±{uj) by adopting some specific driving 
terms g{X) appearing in the main text. First let us determine the factors G±{uj). Using the 
explicit form of K{(jj) in ([H]), we have 

2e-"l'^l/2posh(ca;/2) = G+{uj)G^{uj). (Bll) 
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The well-known relation r(l/2 + a;)r(l/2 — = vr/ cosvrx 58| together with the asymptotic 



form r(x) ~ y2nx^ ae ^ for |x| ^ 1 and the condition for G±{ijj) in Eq. (jBSj) yield 

2n ( 2716 \ 2- 



G+iuj) = G(-io) = — — — (B12) 



where e — )■ 0+. Useful special values are 











(-?) 


T 







G,(0) = v/2, G,(-)=G_(--)=J^. (B13) 



1. The case g{X) = o + bs{\ + Aq) 

Set the driving term to be g{\) = a + hs{\ + Aq), where a, 6 G C. Taking the Fourier 
transform yields 

g{u) = 27ca6{u) + — — -. (B14) 

2cosh(cco'/2) 

Let us decompose the function as in (IBSp . The first term can be easily decomposed by using 

S(uj) = ^(^. —] (e^+0). (B15) 



27ri \uj — ie uj + ie ^ 

On the other hand, the second term in (1B14P is a meromorphic function whose poles are 



simple poles (denoted by ujn) located at 



C 

Thus the decomposition of the function 1/ cosh(cw/2) reads 

1 



a;„ = — (2n + l) (n G Z). (B16) 



cosh(caj/2) 



oo 



2^ ^ (-1 



cosh(ca;/2) c a; + a;„ 



Using this, we can express the function f_{uj)/ cosh(ca;/2), where f^{uj) is any function 
which is analytic and bounded in the lower half-plane, as 

TT 77^ = 

cosh(ca;/2j 

2z ^ 



c ^ w + 

n=0 



= cosh(cc./2) - 7 ^ — ^tt:;:^ — • 

31 



Applying the formula (IBISP to (1B14P and (IBSp . we obtain, for instance, 



n=0 



Substitution of the above equation and ( 1B16I) into ( 1B9I) then yields 



a ^— ^ + b ^ + 

00 + ze c 00 + m/c 



(B20) 



2. The case ^(A) = aK{X) 



Next we consider the case g{X) = aK{\) (a G C). Using the factorization equation ( IBSp . 
one finds 



Thus we immediately obtain defined by Eq. (IBSp . namely 



*+M = -7rTT' $-M = «G_M. (B22) 



y+{uj) = a{G+{uj) - 1), y_{uj) = a ( 1 - -i— ) . (B23) 



Then ^±(0;) are respectively given by 

Appendix C: Low-temperature behavior for arbitrary c > 

In Section III, we derived the explicit low-temperature thermodynamics (see (H2|) for in- 
stance) for the strong coupling regime (c ^ 1) in weak magnetic field {H ^ 1) via the 
polylogarithm expression for the pressure (1371) . Here, using the dressed function formalism, 
we extract the universal low-temperature thermodynamics (jinj) for arbitrary repulsive cou- 
pling c > in arbitrary magnetic field H < He- As shown in Section [XTTt the low-temperature 
thermodynamics is characterized by the two integral equations 

e{k) = k^-^-^-Tai* ln(l + e-'^'^^^^^) 
0i(A) = H-Tai* ln(l + e'^^^)/^) + Taa * ln(l + e-^^^^)/^). (CI) 

In the limit T — )■ 0, the above equations coincide with the dressed energies (lA4p . 
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In completely the same way as in the derivation for Eqs. (fT6|) . ( 123|) and (12^ one obtains 
£(k) = e^^\k) + / ° ai{k - A)0i(A)rfA 
^^(X) = ^f>+ ai{\ - k)e{k)dk - a2{\ - ^i)<pMd^ (C2) 



{ai{k - Ao) + ai{k + Aq)) 



where 








= k' 


H 




= H 





(ai(A - A;o) + ai(A + fco)) + ^ . MX - Aq) + a2(A + Aq)) . (C3) 



6£i(Ao 

In this limit, the free energy (|9]) is reduced to 

^2 1 /"fco 



ttT^ 1 



(C4) 



-fco 

Applying the dressed function formalism to (lC2p and (]A1|) . we arrive at 



1 

2^ 



e{k)dk= \ e^^^p,{k)dk+ / (t)fps{\)d\ 

fco J —ko J -\o 

This yields 

+ (C6) 

6 \Vs vj 

where Vc and Vg are, respectively, the holon and spinon excitation velocities 

^ 4(fco) ^ <(^o) /r.7^ 

Before closing this Appendix, we reproduce the low-temperature thermodynamics ( 142|1 
for c ^ 1 and H = 0. At H = 0, the Fermi point Aq = oo. By Fourier transformation, the 
density functions Pc{k) and Ps{X) defined by Eq. (lAip reduce to 



"fco 

Pc{f^) = 7^+ I ai{k - X)ps{X)dX 

-fco 



1 

2^ 

/fco 
s(A - k)p,{k)dk. (C8) 
■fco 



Up to leading order in 1/c the second equation reads 



ffco 

p,(A)^s(A) / p,{k)dk = n,s{X). (C9) 

-fco 
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Substituting this equation into the first equation in Eq. (ICSP yields 

1 f°° 1 1 In 2 
Pc{k) ^— + nj ai{k- X)s{X)dX = — + n,K{k) ^ — + (CIO) 

where ( 135|) has been used. The dressed energies can be obtained by just taking the hmit 
T ^ 0, ^ and Ao ^ oo in Eqs ([36]) and (Il6]). Exphcitly, 



e^k) ^k"^ - II- 27rPoK{k) ^ k"^ - ii - 2Po In 2/c 

SsiX) ^ -2nPos{X). (Cll) 



The pressure Pq is determined by solving f lA6p . Combining Eq. fIClOP with Eq. flA2p gives 
the Fermi point 

A 21n2\ 

fco ^ ^^-^ ^ t™, (^1 - -— j . (C12) 

fco is derived from the condition (lA5p . with result 

9 2Poln2 
A;2 = /i + ^ . C13 

c 

Substituting these expressions into (1A6P and taking terms of order 0(l/c), one arrives at 

2n37r2 / 61n2\ .^^ ^, 

Po = -Eo + liric = (1 - — p J (C14) 

where fi denotes the zero-temperature chemical potential determined by Eqs. (IC13P and 

5 , / 161n2\ , 

/X = vrX (^1 - j • (C15) 

The excitation velocities 

/ 41n2\ 27rV / 61n2\ , 

are evaluated by substituting all the results into Eq. (IC7p . Thus from Eq. (IC6p . one finds 
that the low-temperature free energy F agrees with Eq. ( H2|) at H = 0. 
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Appendix D: Correlation functions for H <^ 1 



Here we consider the leading terms of the zero temperature correlation functions listed 
in Section V. The conformal dimensions are 



and 



+ ± l^N^ - ± (I) D. (d, + \d. ± Iak 



41n2 

7 



4 fH \ [AN, 



(Dl) 



2At{AN,,s.Nl.Dc,s)^ 

\ (aNs - Ian, ±d,] + ^ (^] AN, (an, - Ian, ± 



Un(Ho/H) 



7r2 \K 



Di -[AN,- -AN, 



± 



H 
7 



D,+ -D, 



AN,--AN,±D, 



(D2) 



(i) The leading orders for the field correlator G^{x,t) come from the set of quantum 
numbers {D„Ds) = (1/2,-1/2) and (1/2,1/2). The conformal dimensions corresponding 
to these sets of quantum numbers are 

H\ \n2 f H 



2A+ 



2A: 



2A+ ^ 
2 

2A; = 



9 31n2 3 
16 ~ 47 ^ 
1 31n2 
16 ~ 
1 



H\ 3 
27r2 V^y ^ 2^ V^c 



47 27r2 Vi^, 
2 1 /i/ 



1 fH 



+ 



71^7 yifc 
3 In 2 f H 



7 V^c 



27 \H,) ■ 71^7 V^^c 
41n2 r H 



+ 



71^7 \ H, 



(D3) 
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for {D„ D,) = (1/2, -1/2) and 



2A+ = 

c 


25 
16 


51n2 

+ 

47 


2A; = 


1 

16 


51n2 

+ —A 


2A+ = 







2a: = 


1 

2 " 


2 / H 



H\ 5 f H\ \n2 f H 
+ — 



27r2 V/Zc/ 27 \HJ 77^7 Vi/c 

1 I f H\ 51n2 /// 



27r2 V i/J 27 \HJ 7727 V //< 



7 V-f^c/ 71^7 
for {D^, D,)^{l/2, 1/2). Hence we obtain 

where the exponents are given by 

5 31n2 1 f H\ 2fH\ 2hi2 f H 
Oci ^ - - + — I — I + - I — 1 + 



27 \HJ 7 \HJ TT^-f \H, 

13 51n2 3 fH\ 2 f H\ 6\n2 f H 



27 7r2 \HJ 7 \HJ 7r27 \H, 



1 2 f H\ 3 / H\ 41n2 f H \ , , 

(ii) The leading terms for the field correlator G^{x,t) come from the set of quantum 
numbers {DcDg) — (0, 1/2) and (1, —1/2). The corresponding conformal dimensions are 
2A+ - — - 3 _ _3_ f H \ _ 3_ f H_\ ln2 f H_ 

1 31n2 1 f H\ 1 /if\ 31n2 /// 
2A7 = — — + 



16 47 27r2 \HJ 27 \HJ 77^7 V^c 

1 2 f H\ I f H\ 41n2 / 
2A+ = - + — — + 



' 2 7r2 Vi/c/ 7 \HcJ 71^7 V^^c 
2A; = (D7) 



for {D^,Ds) = (0,1/2) and 



2A+ = - + ^ + — (—\ - - (—] + 1^ 2 / i7 



16 47 2^2 ^ff j 27 V^c/ 7^27 \H, 

1 51n2 1 1 / H\ 51n2 / 

2A- = — + ^ + __+__ + 



16 47 27r2 \HcJ 27 V^c/ 7r27 \H, 
2At = 



1 2 f H\ 3 / H\ 41n2 / \ , , 
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for (Dc, £>s) = (1, -1/2). This gives 



where 



cl 



5 31n2 1 






( 


8 27 7r2 ' 






7r^7 \ 



13 51n2 3 /H\ 2 / \ Q\n2 f H 
1 2 f H\ I f H\ 41n2 / 



^c2 = - + -^ + -( — |--( — 1 + 
^.1 = ^ + ^ ( ^ 1 + 



1 2 3 41n2 , 

(iii) The leading terms for the charge density correlator come from the set of quantum 
numbers {Dc,Ds) — (0,0), (0,1), (1,0) and (1,-1). The corresponding conformal dimen- 
sions are 

2A+ = 0, 2A+ = 0, 

2A;=0, 2A; = 0, (Dll) 



for (De, Ds) = (0, 0) and 



2A+ = - + 



In 2 




{- 


7 






In 2 




{- 


7 







41n2 f H 



1 ln2 2 ( H\ 41n2 / 
2Ar = - + 



2A+ = - -I I 

' 2 ^ AHHo/H) ^ 7 \H, 

1 1 2 f H\ 



for (Dc, D,) = (0, 1) with 



* 4- 4 In 2 , , 

2A+ = 1 + — -, 2A+ = 0, 

4 In 2 

2 A; = 1 + — — , 2 A; = 0. (D13) 
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for (Dc, £>s) = (1, 0) and 



, , 1 ln2 2 f H\ 41n2 / if 
2A+ = - + + — _ + 



1 In 2 2 f H\ A\n2 f H 
2A- = - + + _ _ + 



4 7 TT^ \HcJ 7r^7 \iic 



1 1 2 /if 

2At = - + 



2 Un{Ho/H) 7 V^^c 



for {Dc,Ds) — (1,-1). The correlation function is then given by 



2 Aicofi{2kF\.x) ^2 cos(2/o>|x) ^3 cos(2(A;ir^ + A;i?i-)a;) 



where 



(D15) 



1 21n2 4 /ii\ 81n2 f H 
Oci = - + 



2 7 7r2 Vi/c/ 7^27 \H, 

1 21n2 4: f H\ 81n2 /if 
^c2 = - + + — I — 1 + 



0c3 = 2 + 



2 7 7r2 Viic/ 7r27 V^^, 
81n2 



7 



1 4 
- ' + 21n(iio/ii) + 7Uc, 



2ln{Ho/H) ^\H, 

(iv) The longitudinal spin-spin correlator is similar to the charge density correlator ob- 
tained above. The only difference is that the leading term is replaced by (m^)^. 

(v) The leading terms of the transverse spin-spin correlator are obtained from the quan- 
tum numbers {Dc,Dg) = (1/2,0) and (1/2,-1). The corresponding conformal dimensions 
are 

, 1 ln2 2 f H 

2A+ = - H — 

" 4 1 l\H, 



1 In 2 2 / H 

4 7 'y \H, 



2A+ 



c 

1 1 2 /ii 



2 41n(iio/ii) 7 V^c 



1 1 2 r H\ 
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for {Dc,D,) = (1/2,0) and 



2A+ = 0, 2A+ = 0, 

2A; = 0, 2A; = 2 (D18) 

for {Dc,Ds) — (1/2,-1). The correlation function is 

^ A,cos{{kn + k,,)x) ^ A,cos{{k,,- k,,)x) 
\x + ivct\^^^\x + ivstl"'^ \x + ivstl^"'^ 

where 

1 21n2 



2 7 
- 1 1 I 

= 2. (D20) 



(vi) Lastly we consider the pair correlator. The leading order terms are contributed from 
the quantum numbers (DcDg) = (1/2,0) and (1/2,-1). The corresponding conformal 
dimensions are 

2A; = l-ii^, 2A,-=0, (D21) 

for {Dc,D,) = (1/2,0) and 



, , 41n2 4 f H\ 81n2 f H 
2 At = 1 + ^ 



41n2 4 81n2 / 

2A- = 1 ^ + 



7 TT^ \HcJ 7r^7 V-^< 



2A- = i-±f^V ^ ^81n2|^£ 

1 4. f H\ 1 81n2 /i/ . , , 

for {Dc,Ds) — (1/2,-1). The correlation function is then given by 



X 



Q AcOs{{kF^ + kF^)x) ^ A2 COs{{kFt-kFl)^j 
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where 

5 61n2 



2 7 
81n2 

9c2 — 2 

7 

We have obtained the charge and spin velocities in Eq. ( ]C16P for this regime. An extension 
to the finite temperature correlation functions is straightforward. 



Appendix E: Correlation functions for H ^ He 



The conformal dimensions in this case are 



2A±(AAr^,„Ar±,D,„ 

2 



\/ 1 - ^DANs + 2Nt (El) 



and 



2A±(AAr,,„Ar±,D,,,) = 



7 ( vr^ 



AA'c - - A A", 



AN, - -AN, 



+ 2N 



± 



(i) G^{x,t): The conformal dimensions are 

2A+ 



2 / if 8 / H 



2A~ = 



2A^ 



2A: 



1 2 1 H A H 

4 ~ 7 ^ 2^y~ If, ^ Ife 



4 7 27r 



He 7r7 



(E2) 



(E3) 
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for {D„ D,) = (1/2, -1/2) and 



2A+ = l + Vl-f + A 

ttV He n-f\ H, 



c 










1 

4 + 


2 


3 


7 


2^ 


1 

4 + 


2 

- + 

7 


1 

2^ 



' - - , HAH 



c 



for (Dc, Dg) — (1/2, 1/2). The correlation function is 

Ai cos(2A:ir-^3;) ^ cos((A;Ft + 2A:j.|)a:) 



with 



Oci = 1--W1- — + — Jl- — 
TT V ifc vi"7 V i^, 



c 



Oc2 = I + -WI- — + — — 

TT V He TT-fM He 



1 4 1 H 16 H 

- 2~7~^V ^~^ ^^V^~^ 

1 4 1 / ^ 

2 7 TT V ^^c 



c 



(ii) Gj.(x,i): The conformal dimensions are 



c 



A- 1 2 1 / 4 / 



c 



C 



2A; = 

for {De,Ds) = (0,1/2) and 



■= 4 7 irV n',\l H, 



c 



12 1/ H 4 / 

4 7 TtV He TT-f]l H, 

2A+ = 



c 



" - , H 24 H 

2a: = i /i_ +_ /i_ 

7 TT V 777 V 



c 
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for {Dc,Ds) — (1,-1/2). The correlation function is 

Aicos{2kFix) A2Cos{{2kF^ + kFi)x) 

Lj\(x,t) !~ -. : — j 77 h 77 — j ] 77 — 

\x + iVcW'^'^lx + iVgtn^ \x + iVctr^'^lx + iVstr^'^ 



where 



6r^ — — 

"^ 2 7 

5 4 4 r w 32 r w 

2 7 TT V He TT-fy He 

2 I W 8 I W 

TT V He n-f]/ He 



8 2/ H 24 H 

7 TT V -f/c 7r7 V -He 

(iii) Gnn^ The conformal dimensions are 

2A+ = 0, 2A+ = 0, 



2A;=0, 2A; = 



for (Dc, D,) = (0, 0) and 



2 / H 16 H 
2A+ = 0, 2A+ = l--./l- — + — Wl-— , 

TT V He TT-fM H, 



c 



ttV He 7r'y\ H, 



c 



for (Dc, ^s) = (0, 1) with 



2A+ = 1 + — Jl-— , 2A+ = 0, 
7r7 V -f^c 

2A; = i + — ./l-|^, 2A; = 
7r7 V -f^c 



for {De,Ds) = (1,0) and 



. 4 / if 16 / H 
2A+ = 1 - -Wl - TT + — V 1 - TT 

TtV He 7T-f\ H, 



c 



4 / H 16 H 
2A~ = 1 A/l ^ Wl 

TtV He^TT^y H^ 



c 



8 2/ H 32 H 
2A+ = 1 Wl \ Wl 

' 7 ttV He n-fy H, 



c 



8 2/ H 32 H 

2A~ = 1 Wl ^ Wl 

' 7 ttV He n-f\ H, 



c 
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for {Dc,Ds) — (1,-1). The correlation function is then given by 

2 , Aicos{2kFix) A2Cos{2kF^x) A3Cos{2{kFi + kF^)x) 

Ljjijiix, t ) ~ 71 + — j r "i r2 — j r i I il/lo ) 

\X + IVst]^"'^ \X + IVct]^''^ \X + lVst\^'^ \X + lVct\^^^ 

where 



= 2-^\/i-#+-\A-# 

TT V i/c 7r7 V H, 

7r7 V 



c 



c 



4 / H 32 H 
7^ V i/c 7r7 V i^c 



es2^2 Ji - — + —Ji - —. EI6 

7 TT V i^c 7r7 V i^c 

(iv) G^{x,t) has the same form as Gnn{x,t) except the term is replaced by (m^)^. 

(v) G'-'-(a;, i): The conformal dimensions are 



1 2 4 
2A+ = - - - + — 

4 7 7r7 

12 4 

2A; = 7 + - + — 
4 7 7r7 

.+ 121 





H 


1 - 




/ 








1 


1 






/ 


1 







for (Dc^s) = (1/2,0) and 

2A+ 



12 2 

4 7 TT 

12 2 
2A; = - + 

4 7 TT 





12 

7r7 \ 



















2A+ 



12 3/ 1 16 / 
4 ~ 7 ~ 2^V ^ ~ ll,^ 



c 



9 6 3 / 1 32 / 
2A7 = 4/1 \ 4/1 (E18) 



for {Dc,Ds) — (1/2,-1). The correlation function is 

G\x, t) ^ + + (E19) 

\x + iVctr'^'^lx + iVstl^^^ \x + iVctr^^lx + iVstr""^ 
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where 

1 S I H 



del = - + — Wl- „ 
2 7r7 V He 

1 4 I W 8 / F 

2 TT V He TTJ \ H, 

1 1 / 



c 



5 8 3 I H 48 I H 

Os2 = ^ — + — Wl- — 
2 7 TT V He 7TJ \ ii, 



c 



(vi) Gp{x,t): The conformal dimensions are 



2A+ 



2A; 
2At 



9 
4 


6 

7 


12 
+ — 

7r7 
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2 


4 


4 


7 


7r7 
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2 

+ - 
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~ 2^ 
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2 




4 


7 


~ 2^ 



for {De,Ds) = (1/2,0) and 



/ 










H 




K 


1 


H 




H'e 


/ 


H 




H'e 



H 

7r7 V H( 



c 



2A: = ^-^-^,/i-4^ + -yi-| 



« 4 7 jrV /fe tttV -ff, 



C 



OA- 1 2 2 r iy 12 r 



c 



1 2 1 / 8 / H 

'^^ = i-^+2^r-H^^r-H: 

OA- - ---- - [7^^ — ll - — 

' ~4 7 2^V K^^y 'H'e 

for (I?c, -Ds) = (1/2, —1)- The correlation function is then given by 

AiC0s((fci7-| + kFt)x A2 cos((/cFt - kFi)x) 



\x + ivetl^^^lx + iVs^l^si \x + ivetl^^^lx + iVst\^''^ 



where 
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^.2 = 7^ W1- — + — W1-— . 

2 7 TT V He nj \ if. 



c 
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Finally we note that the charge and spin velocities can be derived easily from the relations 

e'iko) <f>[{Xo) 



The leading terms in the velocities are then found to be 



(E26) 



v,^27rnAl-—Jl-^] (E27) 



7r7 V H^ 



c 



and 
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